Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



O 

o 

(N 

>>: 



(N 



I 

o 



00 

cn 
in 
o 

o 



I 

O 

o 



X 
J3 



Spin- Wave Theory of the Multiple-Spin Exchange Model on a Triangular Lattice in a 

Magnetic Field : 3-Sublattice Structures 

Chitoshi Yasuda*, Daisuke Kinouchi and Kenn Kubo 

Department of Physics and Mathematics, Aoyama Gakuin University, Sagamihara 229-8558, Japan 

We study the spin wave in the S = 1/2 multiple-spin exchange model on a triangular lattice 
in a magnetic field within the linear spin-wave theory. We take only two-, three- and four- 
spin exchange interactions into account and restrict ourselves to the region where a coplanar 
three-sublattice state is the mean-field ground state. We found that the Y-shape ground state 
survives quantum fiuctuations and the phase transition to a phase with a 6-sublattice structure 
occurs with softening of the spin wave. We estimated the quantum corrections to the ground 
state sublattice magnetizations due to zero-point spin-wave fiuctuations. 

KEYWORDS: multiple-spin exchange model, three-sublattice structure, spin-wave approximation, soften- 
ing, sublattice magnetization, quantum phase transition, solid helium, frustration 



1. Introduction 

The multiple-spin exchange (MSE) model on the tri- 
angular lattice is one of the two dimensional frustrated 
quantum spin systems which are presently at the fo- 
cus of strong interest. The model is believed to describe 
nuclear magnetism of the solid '^He layers adsorbed on 
graphite, ^^'^ whose peculiar temperature dependence of 
the specific heat with a double-peaked structure is quite 
intriguing.^ Several studies were devoted for theoretical 
understanding of the experimental results in terms of 
the MSE model. ^"^'^ The mean-field theory was applied 
to the MSE model with two-, three- and four-spin cyclic 
exchange interactions, which lead to various kinds of the 
ground state corresponding to different values of parame- 
ters.^' Appearance of various mean-field ground states 
reflects weak stability of the mean-field ground states 
due to strong frustration. Therefore quantum fluctua- 
tions should have strong influence on the ground state 
of the system. Quantum mechanical ground state was 
studied by the numerical diagonalization of finite clus- 
ters and a spin liquid ground state with a finite spin gap 
was predicted for the parameters which are thought to be 
adequate to the '^He layer. However recent measure- 
ment of the susceptibility down to ~10 jiK did not show 
any sign which suggests the existence of a spin gap.^"' At 
present, the ground state of both the two-dimensional 
solid '^He on graphite and the MSE model on the trian- 
gular lattice are, therefore, hardly understood. The MSE 
interactions might be relevant as well to the spin liquid 
state of K-(ET)2Cu2(CN)3 where antiferromagnetic in- 
teraction is thought to be dominant. Therefore it is 
meaningful and desired to study the properties of the 
MSE model on the triangular lattice into detail. 

In the present work, we study the spin wave in the 
S — 1/2 MSE model with two-, three- and four-spin ex- 
change interactions on the triangular lattice in the mag- 
netic field. The system is described by the Hamiltonian 



where cTj is the Pauli matrix on the site i and the sum- 
mations and run over all pairs and minimum 
diamonds, respectively. Since the three-spin exchange in- 
teraction is reduced to the conventional two-spin one, the 
first term of the Hamiltonian (1) describes the sum of 
two- and three-spin exchange interactions. The Hamilto- 
nian hp is the four-spin exchange interaction for a min- 
imum diamond p. For a diamond of spins 1 4 with 
diagonal bonds (1,3) and (2,4), it reads 



l<i<j<4 



(Ti ■ cr-j + (cTi • cr2)(cr3 • 0-4) 



(1) 



<»,i> 



+ (cTi ■ cri){a-2 ■ 0-3) - (cTi ■ o-3)(cr2 • 0-4) . (2) 

The coupling parameters are written in terms of con- 
ventional exchange constants as K — — J4/4 and J = 
— J2I2. Since J„ are known to be negative in solid 
^He,"'^^ we take the value of K positive in the following. 
The third term in Eq. (1) is the Zeeman term with a 
magnetic field h [h = \h\) applied anti-parallel to the 
z-direction as shown in Fig. 1. 

Numerous mean-field ground-state phases were found 
for the Hamiltonian (1) by assuming up to 144 sublat- 
tices.^^ The phase diagram is parametrized by J / K and 
h/K. When ft, = 0, the ground state varies according to 
J/K as follows : I) the ferromagnetic phase for J/K < 
—8.61. II) the intermediate phase for —8.61 < J/K < 
—2.26 where the ground state spin structure varies with 
the change of J/K and ten small phases with equal to 
or more than 12 sublattices were identified. These small 
phases might be artifacts of assumptions of finite num- 
ber of sublattices. The true ground state in this param- 
eter region is still controversial even at the mean-field 
level. Ill) the tetrahedral phase for -2.26 < J/K < 8.22. 
The ground state has a 4-sublattice spin structure with 
zero magnetization, where spin vectors on four sublat- 
tices point four vertices of a tetrahedron if their bottoms 
are put at its center. IV) a phase with six sublattices for 
8.22 < J/K < 10. The groimd state has non-coplaner 
spin structure with uniform vector chiral order and stag- 
gered scalar chiral order. Finally, V) the 120° phase for 
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Fig. 1. Y-shape state on a triangular lattice in a magnetic field 
h. 



J > lOK. The spins on three sublattices make the angle 
2tt/3 to each other. 

The mean-field phase diagram in the magnetic field 
was also studied.^' ^^■^'^ One of the interesting behaviors 
is the appearance of a magnetization plateau with 1/2 
of the full polarization in the parameter region adequate 
to '^He layers.^^^'^^ This plateau is realized by the 4- 
sublattice uuud spin structure. Other numerous ground- 
state phases in the magnetic field have been found, but 
we introduce only the phases related to the present work. 
The 120° ground state for J > lOA' is modified by the 
magnetic field. This state has a 3-sublattice coplanar spin 
configuration with the spins on one (say A) sublattice 
antiparallel to the applied magnetic field. The spins on 
other sublattices (B and C) tilt toward the oblique direc- 
tions as shown in Fig. 1. In the present paper, we call it 
Y-shape state. This result shows that the four-spin ex- 
change interaction lifts the non-trivial degeneracy of the 
mean-field ground-state of the antiferromagnetic Heisen- 
berg (AFH) model in the magnetic field. A 6-sublattice 
phase exists for 8 < J/K < 10 adjacent to the Y-shape 
state for weak magnetic field. Increase in the magnetic 
field induces a 12-sublattice structure. Further increase 
in the magnetic field leads to a magnetization plateau 
with 1/3 of full polarization due to the 3-sublattice uud 
state with two sublattices with up spins and the other 
with down spins. The uud state is stabilized by the four- 
spin interaction and occupies a large region in the phase 
diagram. We note that the phase boundary between the 
Y-shape and the 6-sublattice phases was not determined 
by the mean-field theory. 

Spin waves in this model were studied previously in 
the absence of the mag netic field.<'''^'i2 xhe spin wave m 
the tetrahedral phase was studied and the quantum cor- 
rections to the sublattice magnetization were estimated.® 
The analysis for J = showed the stability of the tetra- 
hedral state against the zero-point fluctuations of the 
spin wave. The spin wave in the 120° phase was studied 
and it was shown that there are three gapless branches 
for small k and the spin wave softens at J = lOKJ The 
softening corresponds to the phase transition to the 6- 
sublattice phase. 

In the present work, we investigate the spin wave in 
the Y-shape phase and discuss the phase transition under 
the magnetic field. Although the Y-shape phase does not 
correspond to the model for the solid ^He layer, coplay 



of the geometrical frustration and the MSE interactions 
poses a problem of interest in itself. We investigated the 
dispersions of the spin waves as well as the effects of the 
zero-point fluctuations of the spin wave on the ground- 
state energies. As a result, we found that though the 
Y-shape ground state survives quantum fluctuations for 
J/K > 12 and small h/K, softening of the spin wave 
leads to the phase transition to the 6-sublattice phase. 

The present paper is organized as follows. In § 2, we 
present the linear spin-wave theory for the Y-shape state 
in the magnetic field. In § 3, we simply summarize the 
magnetic properties in two special cases, i.e., the AFH 
model in the magnetic field and the MSE model at /i = 0. 
The dispersion of the spin wave and the phase diagram 
obtained from the softening of the spin wave are shown in 
§ 4 and 5, respectively. The quantum effects are discussed 
by studying quantum corrections of the ground-state en- 
ergy and the sublattice magnetizations in § 6. Finally, 
§ 7 is devoted to summary and discussion. 

2. Spin- Wave Theory for the Y-Shape State 

Assuming the Y-shape state as the ground state, we 
can perform the Holstein-Primakoff transformation^*^ of 
the Hamiltonian (1) by neglecting higher-order terms as 




(3) 



for i e A sublattice. 



a^^-a{bl + b,)-(3{l-2b%) 
at ~ ijafl - 2b]b,) -{13 + l)b] -{(3- 1)6,} 



l{a{l~2b%)-{(3-l)b]-iP + l)b,} 



(4) 



for j G B sublattice, and 

<jlo,a{ci+Ck)-(3{l-2clck) 
a+ i{-a(l - 24c,) - (/3 + 1)4 - (/3 - l)c4 , 
a- c i{-a(l - 24cfe) -if3- 1)4 - (/3 -f l)c4 

(5) 

for fc e C sublattice, where a- (a^), 6| {bj) and 4 (cfc) 
are the boson creation (annihilation) operators, a — sinO 
and /3 = cos 6 as shown in Fig. 1. 

First, we obtain the relation between h and the angle 
9 from the condition that the first-order terms of the 
transformed Hamiltonian should vanish, i.e., 

1 



12K 
for h < he where 



{5K + J- ^/{J - Kf - AKh} , (6) 



he = -12K + 3J. 



(7) 



The value of monotonously decreases with the increase 
of h and the uud state {6 = 0) is realized for h > h^. This 
relation can be also obtained from dE'^^/dO — where the 
mean-field ground-state energy iJ'^' per the total number 
of sites N is given by 

E^^/N = {2(3^ - 2/3 - 1) J (8) 
-(8/3^ - 10/32 + l)K + hi - 2p)h . 
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Fig. 2. Unit vectors <5i = (1, 0) and 62 = (0, 1) on the triangular 
lattice. 



After straightforward calculations the Hamiltonian is 
rewritten as 

^ = ^" + E^fe^^^fe+^o , (9) 

k 

I 

where v^"^ = (a^, 6^, c^, a_fc, c_fc) and ^ in the 

k 

second term denotes the summation over a half of the 
reduced Brillouin zone of the 3-sublattice structure. The 
matrix V reads 



V 



Moff Mdiag 



where 



and 



with 



Mdiag = ( Cr* Bk DTk 
CTk DTI Bk 



M, 



off 



^ GTk HTl, Fk 



(10) 



(11) 



(12) 



Ak = 12{-2K +{J + 2K)(3 + 4^X13^ 

+K{1- (3^)Ak]-2h , 
Bk = 6{ J + 2K+{ J - 4:K)(3 - 2{ J + 5K)f 

+ 12K/3^ + 2K(1 -fi + 2f3^)Ak} + 2h(3 , 
G = 3(1 - /3)( J + AK + AK[3^) , 
D = 6{X(4/3 - 1) + ( J + bK)0^ - 8(3^ K} , 
Ek = -12{l~p^)KAk , 
Fk = 2f3Ek , 

G = -3(1 + f3){J + 4f3^K) , 
H = -6(1 - f3^)iJ + 5K~ 8PK) . 
Here the wave-number dependence is incorporated by 



(14) 



Afc = 3't'^°^ ('^i + /!;2)+cos (2fci - fc2)+cos (-/ci + 2fc2)} , 

(15) 

where fc^ = fe • 5^ j = 1 , 2 is an element of the wave vector 
and the unit vector and 62 of the triangular lattice 




Fig. 3. Three dispersions of the spin wave for the AFH model 
with (a) /i = and (b) h/J = 2. The marks F, M and K denote 
the wave numbers as defined in Fig. 4 and the vertical axis the 
Unes F-M-K-r. 



are chosen as shown in Fig. 2. Note that kf - 
3fcV4 holds. The third term in Eq. (9) 



F'^ 



k 



(Ak + 2Bk) 



kik2 = 



(16) 



arises from the Bose commutation relation. The Hamil- 
tonian (9) is diagonalized through the Bogoliubov trans- 
formation and takes the following form: 



F'' + Y,Uk^Vdu 

k 



,Uk 



(17) 



where Uk^ = (aj., 7]., Q;_fc, 7-fc) is a vector of 
the transformed Bose operators and 2?diag is a diagonal 
matrix whose elements are lo^^ 



(3) 

(13) and ujj^ for « = 1 ^ 6. The frequencies of the three 

branches of the spin wave are given by w^^"* , uj^^"^ and w^^'' . 
We numerically performed the transformation according 
to the general theory by Colpa.^^ 



, ,(2) , ,(3) 
''^fc ^^k ' 



(2) 
k 



Spin- Wave Spectrum of the AFH Model in 
the Magnetic Field and the MSE Model at 



In this section we reproduce the previously known re- 
sults for simple cases in order for comparison with the 
results in the next section. 

First we demonstrate the results for K = 0, i.e., the 
AFH model. The dispersions of the spin waves assum- 
ing the Y-shape state as the ground state for the AFH 
model with h/J = and 2 are shown in Figs. 3 (a) 
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Fig. 4. Brillouin zone of the 3-sublattice structure with reciprocal 
lattice vectors (fci, fc2) = (27r/3, -27r/3) and (27r/3, 47r/3). The 
marks P, M and K denote (fci, ^2) = (0,0), (27r/3, 7r/3) and 
(27r/3, 27r/3), respectively. 



and (b), respectively. The hexagonal Brillouin zone for 
the 3-sublattice structure is shown in Fig. 4. For h = 
the frequencies of all spin-wave branches vanish at the F 
point (fc = 0) and have the linear dispersion for small k. 
In the magnetic field, one of the three branches is lifted 
to ujk = 2h at the F point corresponding to the uniform 
precession of the spins about the magnetic field. There 
are two gapless branches rather than one expected from 
the SO (2) symmetry. Furthermore one of them has a 
quadratic dispersion for small k. This is caused by a non- 
trivial continuous degeneracy in the mean-field ground 
state in the magnetic field. Any three-sublattice spin 
structure fulfilling the relation Sa -I- S'b -I- S'c = h/3J 
is a ground state for h < 3J, where S'a.b.c is the sub- 
lattice magnetization on each sublattice. Chubukov and 
Golosov showed that quantum fluctuations lift this non- 
trivial degeneracy and select the coplanar spin structure, 
i.e., the Y-shape state. This is an example of the so- 
called order from disorder phenomena. 

Next we show the result for the MSE model at h = 
0. In this case the transformed Hamiltonian (9) can be 
analytically diagonalized^ as 

fc 

(18) 

where the summation of the second term runs over the 
Brillouin zone, and 



4"' = V^i'^ - 



Vl'> = 6( J -K) + 9KAk + 3(J + 5K)r'k 



and r; = Re(Ffc), K 



.(2) 



y(l) Ty(2) 



(19) 

(3) _ 



V-'k+Q, W^f = W'-fc+Q where Q = (-27r/3, 2^/3). The 
first and third terms of Eq. (18) are 



--{J + K)N 



M = l k 



(20) 



(21) 



We show the spin- wave spectrum for J/K = 12 and 10 
in Figs. 5 (a) and (b), respectively. The frequencies of 
branches vanish at the F point as for the AFH model. 




Fig. 5. Three dispersions of the spin wave for (a) J/K = 12 and 
(b) J/K = 10 at /i = 0. The marks P, M and K denote the wave 
number vectors as shown in Fig. 4 and the vertical axis the lines 
P-M-K-P. 



They have linear dispersions for small k as 



,(1) 



,(2) 



3^/3{J - K){J + 2K) k , 

, ,(3) 



(22) 
(23) 



~ 3^6{J - K){{J - lK){kl + kl) - (J + 5K)kik2} . 

Figure 5 shows that the frequency of the second branch 
LOf^ vanishes at the M point, which indicates that a 
second-order phase transition to the 6-sublattice struc- 
ture occurs due to the softening of the spin wave just 
at J = IQKJ The softening occurs at six equivalent M 
points in the Brillouin zone. 

4. Spin- Wave Spectrum of the MSE Model in 
the Magnetic Field 

In Figs. 6 (a) and (b) we show the dispersion for 
J/K = 12 and 10.14 at /i = lOK. In contrast with the 
spectrum in the AFH model, there is only one gapless 
branch at the F point and the gapless branch has a lin- 
ear dispersion. The result is consistent with the excita- 
tions from a stable ground state only with the global 
S0(2) symmetry and corresponds to that the four-spin 
exchanges stabilize the coplanar Y-shape state already 
in the mean-field approximation. The frequency of the 
lowest branch decreases at the M point with decreasing 
J/K for a fixed h. It vanishes at a critical value of J/K 
as shown in Fig. 6 (b). The J/iC-dependence of the fre- 
quency of the spin wave at the point M (fc = (27r/3, 7r/3)) 
is shown in Fig. 7 for various values of h. The softening of 
the spin wave is induced by the competition of the two- 
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Fig. 6. Three dispersions of the spin wave for (a) J/K = 12 and 
(b) J/K = 10.14 at h/K = 10. The frequency of the lowest 
branch at the M point vanishes at J/K C:^ 10.14. 



and four-spin interactions in the magnetic field. The fre- 
quency vanishes apparently as oc (J - Jcf. Here we 
assume the critical exponent /3 to be 1/2 which results 
from the general Landau theory. The softening of the spin 
wave leads to a phase transition to the 6-sublattice phase 
in the magnetic field. We determine the phase boundary 
between the Y-shape and the 6-sublattice phases by de- 
termining the critical value in the next section. 

5. The Phase Diagram 

In Fig. 8 we show the phase diagram determined from 
the softening of the spin wave together with the phase 
boundaries determined by the mean-field calculations. 
The phase boundary between the Y-shape and the 6- 
sublattice phases extends from {J/K, h/K) = (10,0) to 
the multi-critical point ((J/ii')muiti, {h/ K)ynuiti) where 
four phase boundaries appear to merge. Phase transi- 
tions from the Y-shape phase to the 6-sublattice and 
the uud phases are continuous. On the other hand, the 
transition between the 6-sublattice and the 12-sublattice 
phases and that between the 12-sublattice and the uud 
phases are discontinuous. The phase boundary between 
the Y-shape and the uud phases is given by the rela- 
tion (7) and the multi-critical point is estimated to be 
((J/i^)muiti, (/i/if)muiti) ^ (11.37, 22.11). Though in 
the AFH model the uud phase occurs only at a criti- 
cal value he = 3J in the mean-field theory, the phase 
extends to a finite region of h in the MSB model. The 
four-spin exchange stabilizes the uud state. It was shown 
that the thermal or quantum fluctuations stabilize the 
uud state. "'^^'^^ The four-spin exchange has a same effect 
and as a result magnifies the 1/3 magnetization plateau. 
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Fig. 7. J/ZC-dependences of the frequency of the spin wave at the 
point M in the magnetic field h/K = 0, 10, 12, 14, 16, 18, 20 and 
22 from top. The solid marks for a;(27r/3, 7r/3) = denote the 
phase-transition points. The bold line for h/K = is analytically 
obtained by using Eq. (19). The broken lines are obtained by the 
least-squares fitting with the fitting function uj^ {J — Jc)^ with 
13 = 1/2. 




Fig. 8. Phase diagram parametrized by the two-spin interaction 
and the magnetic field. The abbreviations of Y-shape, 6-sl, 12-sl 
and uud denote the phases with the Y-shape, 6-, 12-sublattice 
and uud structures. The open squares are estimated by the soft- 
ening of the spin wave at fc = (27r/3, 7r/3). The diamonds denote 
the parameters at which the total magnetization is equal to the 
saturated value 1/3 of the uud structure. The broken line is the 
phase-transition line h^ = — 12_fS' -|- 3 J. For reference, the data in 
ref. 13 estimated by the vector and scalar chiral order parameter 
within the mean-field approximation are also shown as circles 
and triangles. The solid square denotes the multi-critical point 
((J/i^)multi, ('i/i<')multi) ^ (11.37, 22.11). The bold line is a 
guide for the eye. 



In the phase diagram, we also mark the phase tran- 
sition points among the Y-shape, 6-, 12-sublattice and 
uud phases, which are estimated by the vector and scalar 
chiral order parameters within the mean-field approxima- 
tion. 

Peculiar properties are expected at the multi-critical 
point. Since the spin-wave analysis is broken up just at 
the multi-critical point, we investigate systems close to 
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Fig. 9. Three dispersions of the spin wave for (a) J/K = 12 and 
(b) J/K = 11.37 (~ (J/i<')multi) at h/K = 24. The flat branch 
is observed along the F-M line. 




20 25 30 35 40 45 50 
h/K 

Fig. 10. Magnetic-field dependences of the frequencies of the low- 
est branch at F (open) and M (soUd) points for J/K = 10, 11.37 
(~ (i^/i^)multi) and 12. All the lines are guides for the eye. 



the point. In Fig. 9 (a) we show the dispersion of the 
spin wave for J — 12K and h = 2AK just on the phase- 
transition hne between the Y-shape and the uud phases, 
which is calculated by the spin- wave theory assuming the 
uud state as the ground state. We found that the lowest 
gapless branch has a tendency to be flat along the F-M 
line. Actually, for J/K = 11.37 ~ (J/J^)muiti the branch 
along the F-M line is flat as shown in Fig. 9 (b). As the 
magnetic field is increased from the value on the phase- 
transition line, the frequency of the flat lowest branch 
for J/K = 11.37 ~ {J/K)i-nuiti linearly increases with 
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Fig. 11. Magnetic-field dependences of the classical part of the 
energy E'^ and the total energy E^°^'^^ per bond normalized by 
J + K. For comparison, the data of the AFH model {K = 0) are 
also shown. The lines on the marks are guides for the eye. 
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Fig. 12. Magnetic-field dependences of the quantum correction of 
the energy per bond normalized by J + K. The lines are guides 
for the eye. 



h/K in the range 22.11 < h/K < 34 as demonstrated in 
Fig. 10. Since the second branch at the F point as shown 
in Fig. 9 (b) crosses the lowest branch at h/K ~ 34, the 
/i-dependence of the frequency of the lowest branch at the 
F point has an inflection at h/K ~ 34. When the value of 
J/K is shifted from the multi-critical point, the flatness 
due to the balance of the two- and four-spin exchanges 
is broken up, though the linearity still holds. In the uud 
phase the frequency of the lowest branch has the relation 
uJk = 2{h - hi) and 2(/i2 - h) for hi < h < {hi + /i2)/2 
and {hi + /i2)/2 < h < h2, respectively, where hi and 
/i2 are the magnetic fields of the lower and the upper 
phase boundaries of the uud phase. Fhe linearity of the 
/i-dependence of the frequency has been also obtained in 
the uud phase for the quantum AFH model. 
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6. Ground-State Energy and Sublattice Magne- 
tization 

The total energy of the ground state is evaluated by 
E"^ + E'^, where E'^ is the classical energy of Eq. (8) and 
can be analytically calculated, and E'^ is the quantum 
energy corresponding to, e.g., Eq. (21) for h = 0. The 
magnetic-field dependences of E^ and E'^ per bonds are 
shown in Figs. (11) and (12), respectively, where we nor- 
malize them by using J -I- to compare with the values 
of the AFH model. It is found that, as J/K is decreased, 
the quantum correction of the energy is increased. 

The total magnetization is defined by using the expec- 
tation value of the spin z component, i.e., 

1 



M ■ 



N ^ 

all i 



(24) 



Note that we take the direction of the magnetic field anti- 
parallel to the z direction. The sublattice magnetizations 
of the A and B sublattices are defined by 

^A = -^Y.^<jt) = -l + AmA, (25) 



N 



ieA 



niB 



respectively, where from Eqs. (3) and (4), 



(26) 



(27) 



(28) 



From a symmetry of the present system, the sublattice 
magnetization of the C sublattice is same as that of the 
B sublattice. The total magnetization is written from 
Eqs. (24)-(26), i.e., 

cl 



AM 



where 



2/3) 



2/3AmB) 



(29) 



(30) 



(31) 



The magnetic-field dependences of M for various val- 
ues of J/K are shown in Fig 13. When h — 0, the total 
magnetization is to be zero since the ground state has 
the 120° structure for all systems. For the AFH system, 
the value of — Af is increased with h and saturates to 
—M = 1/3 at h/J — 3, where the phase transition to 
the uud phase occurs. While we can obtain the same 
behavior as the AFH system for J/K = 14, the val- 
ues of —M become negative for small h/{J + K) for 
10.5 < J/K < 12. This result shows that the spin-wave 
theory on the sublattice magnetization is broken up due 
to strong quantum fluctuations in the region. We also 
show the /i-dependencc of the quantum correction AA/ 
in Fig. 14. The quantum corrections vanish both aX h = Q 
and the phase transition point to the uud phase. 




Fig. 13. Dependences of the total magnetization on the magnetic 
field normalized hy J+K. The data described as AFH are results 
for K = 0. 




Fig. 14. Dependences of the quantum correction of the total mag- 
netization on the magnetic field normalized by J + K. 



The magnetic-field dependences of the sublattice mag- 
netizations of the A and B sublattices for various values 
of J/K are shown in Fig. 15. The values of — wa and 
me are decreased with J/K and h/{J + K) except near 
the phase transition point to the uud phase. The reduc- 
tion shows that the quantum effects are enhanced by the 
addition of the four-spin interaction in the MSE model. 
We also show the dependence of the quantum corrections 
of the sublattice magnetizations on the magnetic field in 
Fig. 16. The quantum corrections are larger for smaller 
J/K and larger h/{J + K). We can obtain an interest- 
ing tendency of the size relation between AmA and Ame 
within the present spin- wave approximation. In the AFH 
system with K — 0, the value of Atoa is larger than that 
of Atob except for /i = 0. The result suggests that the 
sublattice magnetization along the axis anti-parallel to 
the direction of the magnetic field is most sensitive to 
quantum fluctuation. In the MSE system with the large 
four-spin interaction and small magnetic field, on the 
other hand, the sublattice magnetization whose quanti- 
zation axis tends to parallel to the direction of the mag- 
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Fig. 15. Dependences of the sublattice magnetizations on the 
magnetic field normaUzed by J + K. The open and sohd marks 
denote — rriA and me, respectively. 



wave to estimate the magnetization in the region. 

The exact diagonaUzation study of finite clusters pre- 
dicted that the ground state is a spin-hquid state with 
a spin gap fiUed with a large number of singlet states 
for parameters corresponding to our 6-sublattice phase. ^"^ 
The character of the ground state in this parameter re- 
gion is, however, still not clarified. In order to numer- 
ically examine the possibility of the 6-sublattice phase, 
we would need the exact diagonalization of clusters with 
larger sizes. The spin-wave analysis assuming the 6- 
sublattice structure as the ground state is another in- 
teresting problem remained for the future work. 
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Fig. 16. Dependences of the quantum corrections of the sublat- 
tice magnetizations on the magnetic field. The open and solid 
marks denote Attia and Arrta, respectively. 



netic field is most sensitive to quantum fluctuation. 

7. Summary and Discussions 

In the present paper, we studied quantum effects on 
the 3-sublattice structures in the S* = 1/2 multiple-spin 
exchange model with two-, three- and four-spin exchange 
interactions on the triangular lattice in the magnetic 
field. By using the linear spin-wave theory, we found that 
the coplanar Y-shape state is stable as the ground state 
of the quantum system, though the four-spin interaction 
leads to the instability as the softening of the spin wave 
and the phase transition to the 6-sublattice phase occurs. 

We had arguments in the framework of the linear spin- 
wave theory. Within the approximation, the magnetiza- 
tion becomes negative in the parameter region where the 
four-spin exchange interaction is dominant. We should 
take into account higher-order corrections of the spin 
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